Transverse momentum broadening of hadrons produced in semi-inclusive 

deep-inelastic scattering on nuclei. 
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The first direct measurement of the dependence on target nuclear mass of the average squared 
transverse momentum {pt) of 7r^,7r~, and K'^ mesons from deep-inelastic lepton scattering is ob- 
tained as a function of several kinematic variables. The data were accumulated at the Hermes 
experiment at Desy, in which the Hera 27.6 GeV lepton beam was scattered off several nuclear gas 
targets. The average squared transverse momentum was clearly observed to increase with atomic 
mass number. The effect increases as a function of and x and remains constant as a function of 
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both the virtual photon energy v and the fractional hadron energy except that it vanishes as z 
approaches unity. 



PACS numbers: 13.60 13.87 14.62 

The evolution of a fast-moving quark into hadrons is 
a non-perturbative and dynamic phenomenon. The ba- 
sic process of hadronization in vacuum is described by a 
well-developed phenomenology, constrained primarily by 
data from inclusive hadron production in ef^ e~ annihila- 
tion. More recently, semi-inclusive hadron multiplicities 
measured in deep-inelastic scattering (DIS) of leptons on 
protons and deuterons, together with inclusive hadron 
yields from the Rhic p-p collider experiments, have pro- 
vided further constraints. For a recent global analysis of 
such data, see Ref. [l[ and references therein. 

The nuclear modification of hadron production in DIS 
was first observed by the pioneering semi-inclusive DIS 
experiments at Slac ^A. followed by measurements per- 
formed by the Emc ll| and the E665 [i| collabora- 
tion. These experiments typically determined the ratio 
of hadron multiplicities observed in the scattering on a 
nucleus to those on Deuterium (D), the so-called nuclear 
attenuation. More recently, much more precise data were 
collected and analyzed by the Hermes collaboration j^- 
as a function of the kinematic variables i^, z, Q^, and 
P( , where v is the energy of the virtual photon, z the 
fractional hadron energy in the target rest frame, — 
being the squared four-momentum of the virtual photon, 
and pt the transverse momentum of the produced hadron. 
It has been suggested that hadron production proceeds 
through three qualitatively distinct stages that involve 
the propagation and interaction of: 

(i) the initial struck quark^ (the "partonic" stage), 

(ii) the subsequently formed colorless state (the "color- 
neutralization" stage, often also termed "pre-hadronic" 
stage), and 

(iii) the final produced hadrons (the "hadronic" stage) . 

This picture is supported by the results of a two- 
dimensional analysis of the multiplicity ratio @ in the 
kinematic variables v and z. However, the existence and 
relative importance of the various stages to the observed 
nuclear attenuation have been difficult to determine un- 
ambiguously. This Letter reports, complimentary to the 
results presented in Ref. the first detailed measure- 
ment of another observable, which may help to better 
constrain models, especially with regard to the role of 
the "partonic" stage: the broadening of the transverse 
momentum distribution of various hadrons. 

In terms of the quark-parton model and QCD, there 
are several contributions to the transverse momentum 
distribution of hadrons produced in semi- inclusive DIS: 



^ "quark" is used in this Letter for both quarks and antiquarks 



(a) primordial transverse momentum, 

(b) gluon radiation of the struck quark, 

(c) the formation and soft multiple interactions of the 
"prc-hadron" , and 

(d) the interaction of the formed hadrons with the sur- 
rounding hadronic medium. 

In semi-inclusive DIS off nuclear targets, the struck 
quark propagates through a "cold" nuclear medium. In 
the nuclear medium the primordial transverse momen- 
tum of quarks may be modified by various effects like 
Fermi motion of nucleons inside the nucleus, modification 
of the nucleon size, formation of non-nucleonic degrees of 
freedom like multiquark states or exchange mesons medi- 
ating the nuclear force. Also, the probabilities of the pro- 
cesses (b-d) may be enhanced resulting in a larger trans- 
verse momentum magnitude of the observed hadrons rel- 
ative to the process in the vacuum or in a free nucleon. 
In particular, process (b) may cause an increased trans- 
verse momentum magnitude and energy loss of the quark 
and therefore it has recently been suggested Q that the 
broadening of hadron distributions in semi-inclusive DIS 
may be the most direct way to probe the "partonic" 
stage. Both the radiative energy loss and the transverse 
momentum of the parton increase with L, which is the 
path length of the quark in the nuclear medium before 
the colorless state formation. QCD predictions for the 
relationship between energy loss and the broadening of 
'Pt distributions are given in Ref. [1, Q. As the rela- 
tionship between these two quantities is independent of 
the dynamics of the initial scattering process, it holds 
equally well in "cold" nuclear matter and in finite-length 
"hot" matter, produced, e.g., in ultra-relativistic heavy- 
ion collisions or high-energy proton-nucleus interactions. 
Thus the understanding of this broadening in the "cold" 
nuclear medium provides precious information for the in- 
terpretation of such high-energy processes [^, • 

An indication of nuclear broadening of the transverse 
momentum distribution for charged hadrons in muon and 
neutrino DIS was reported previously 0, 0, [llj . This 
Letter presents the first direct measurement of the pt- 
broadening for identified mesons and various nuclei as a 
function of the relevant kinematic variables. 

The transverse hadron momentum -pt is defined rela- 
tive to the direction of the virtual photon, which is deter- 
mined by the kinematics of the incident and the scattered 
lepton. The pt-broadening A(pj)^ is defined as the dif- 
ference of the average squared transverse momentum of 
the detected hadron of type h produced on a nuclear tar- 
get with atomic mass number A and that on a deuterium 
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TABLE I: Accumulated yield of selected hadrons for the var- 
ious nuclear targets. 



(D) target: 



MPtA = {Pt)A - (Pt) 
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The broadening may depend on the hadronic variable z, 
and on the leptonic variables v, Q^, and x, with x = ^^j^ 
the Bjorken variable and M the nucleon mass. 

The measurements were performed at Hermes using 
the 27.6 GeV lepton (e+ or e~) beam stored in the 
Hera ring at Desy. The targets consisted of polar- 
ized or unpolarized D, unpolarized He, Ne, Kr, or Xe 
gas. The target atoms were injected into a 40 cm long, 
open-ended tubular storage cell through which the lep- 
ton beam passed. Target areal densities higher than 
10^^ nucleons/cm^ were obtained for unpolarized gases 
and of about 2 ■ 10^^ nucleons/cm^ for polarized D. 
The scattered lepton and the produced hadrons were 
detected in coincidence by the Hermes spectrometer 
Leptons were distinguished from hadrons with an 
average efficiency of 99% and a contamination level of 
less than 1% by using a transition-radiation detector, 
a scintillator pre-shower counter, a dual-radiator ring- 
imaging Cerenkov (Rich) detector, and an electromag- 
netic calorimeter. Scattered leptons were selected by im- 
posing the constraints > 1 GeV^, W'^ > 10 GeV^, 
and y < 0.85, where is the squared invariant mass 
of the photon- nucleon system and y = v/E the energy 
fraction of the virtual photon in the target rest frame, 
with E being the beam energy. In order to suppress the 
contributions from target remnant fragmentation, the re- 
quirement z > 0.2 was imposed. 

The numbers of identified and selected charged pions 
and positive kaons are given in Table [H Negative kaons 
and anti-protons are not considered in this analysis due 
to low statistical precision. Protons were not considered 
because they are not well enough described by the Monte 
Carlo simulation used for the unfolding procedure de- 
scribed below. 

The data used for the analysis of pt-broadening were 
corrected for QED radiative effects, instrumental resolu- 
tion, and acceptance. This was achieved using a Monte 
Carlo simulation based unfolding procedure [l3[, which 
accounted for event migration between different kine- 
matic bins. The Pythia event generator in conjunc- 
tion with a special set of Jetset fragmentation param- 
eters [Til, was used together with the Radgen [T^ sim- 



ulation of QED radiative effects and with a Geant [i3| 
based simulation of the detector. The possible depen- 
dence of the unfolding procedure on the target type 
used in the Monte Carlo simulation was studied with 
the Lepto generator and Jetset using D and Xe 
targets. The results for both targets were found to be 
consistent. Therefore, a Monte Carlo simulation for the 
D target was used to correct the results for all targets. 

Selected events were binned in and in either z, v, 
(5^, or X. These two-dimensional ^ distributions were 
unfolded using smearing matrices that hold information 
about migration among kinematic bins. This method 
inflates the statistical errors to account for this migra- 
tion [l^. The impact of the unfolding on (p^) ranges 
from 60% to almost zero with increasing ly, and is about 
20% for all z bins. It varies from 20% to 30% with in- 
creasing Q^, and is between 15% and 40% for the x-bins. 
After this correction, the average value was calculated 
in each bin of z, x, or and then {p1)d was sub- 
tracted from {p1)a- The charged pion yields were cor- 
rected for pions from the decay of exclusively-produced 
mesons. A correction for such pions, where at least 
one decay pion was in the acceptance, was made using 
a special Monte Carlo generator [2^ that was tuned to 
Hermes data to obtain the correct relative contribution 
from coherent and incoherent p° production This 
Monte Carlo generator also simulated the measured kine- 
matic dependence of coherent production from heavier 
targets. The exclusive p° correction was found to be less 
than 1% over the whole kinematic range, except at the 
highest z-value where it was about 10%. The contribu- 
tion of vector-meson decay to the kaon sample was found 
to be negligible and no correction was included for this. 
A systematic uncertainty due to Rich hadron misidenti- 
fication was estimated to be smaller than 1%. 

The pt-broadening for 7r+ , tt" , and was determined 
for four nuclei as a function of either z, v, x, or , while 
integrating over the other kinematic variables. The aver- 
age kinematics of the p4-broadening results are shown in 
Table im The results reveal negligible correlations among 
the variables, except for x and that are correlated due 
to the forward acceptance of the Hermes spectrometer. 

The nuclear mass dependence of the p4-broadening is 
presented in Fig. [TJ The broadening increases with mass 
number A. It is similar for the charged pions and seems 
to be systematically higher for positively charged kaons. 
The precision of the data does not allow a firm conclu- 
sion about the functional form of the increase of the data 
with A. There is no clear indication of a saturation of the 



^ A two-dimensional unfolding still leaves room for acceptance ef- 
fects if both the yield and the magnitude of the effect depends 
non-linearly on the variables that are not included in the two- 
dimensional binning. An unfolding in more than two dimensions 
was not feasible in this analysis. 
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FIG. 1: The pt-broadening for tt"*", tt~ , and K'^ mesons as 
a function of atomic mass number A. The inner error bars 
represent the statistical uncertainties; the total bars repre- 
sent the total uncertainty, obtained by adding statistical and 
systematic uncertainties in quadrature. 



Pt-broadening at large atomic mass numbers, support- 
ing models which treat its origin in the partonic stage. 
Within such models, this behavior suggests that the color 
neutralization happens near the surface of the nucleus or 
outside for the average kinematics of this measurement 

The panels presented in Fig. [2] show (p^) for D (top 
row) and the pt-broadening (remaining rows) as a func- 
tion of cither Q^, x, and z for tt"*" or 7r~ for the various 
nuclear targets. Since the uncertainties of the sam- 
ple are rather large, only the results for the Xe target 
are presented in the bottom row. The values of (p^) for 
D are between 0.2 and 0.4 GeV^ while the pt-broadening 
shows values from up to 0.05 GeV^. This means that 
Pt-broadening adds between to 10% to (pt). The data 
do not reveal a significant dependence on v in the kine- 
matic range covered. 

Since models that describe hadron formation in nuclei 
commonly connect formation length with z/, the basically 
flat behavior in i/ supports again the picture that color 
neutralization mainly happens at the surface (or outside) 
of the nucleus for the Hermes kinematics [22j- The effect 
slightly increases with in contrast to the model cal- 
culation in Ref. [2^ , where a decrease of the broadening 
with is predicted, and in agreement with the model 
calculation in Ref. 2J]. The behavior as a function of x 
is very similar to the behavior, due to a strong cor- 
relation between x and in the Hermes kinematics, 
hence it can not be excluded that the dependence ob- 
served is actually an underlying x dependcn ce or both 
a and x dependence. The statistical precision of the 
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TABLE II: Average kinematics for the (tt^) pt-broadening 
results. The i/, Q^, and z kinematics are for the Xe target. 



data presented here do not allow the study of the 
and X dependence separately, or any other two kinematic 
obscrvables. 

The pt-broadcning is seen to vanish as z approaches 
unity while the (pj) for D is 0.2 or higher in the high- 
est z-bin. Due to energy conservation the struck quark 
cannot have lost energy when z = 1, leaving no room 
for broadening apart from a possible modification of the 
primordial quark transverse momentum. The observed 
vanishing of the A{pf)\ at high values of z indicates that 
there is no or little dependence of the primordial trans- 
verse momentum on the size of the nucleus. It also indi- 
cates that pf-broadening is not due to elastic scattering 
of pre-hadrons or hadrons already produced within the 
nuclear volume, as this would lead to substantial broad- 
ening even for values of z very close to unity. 

In summary, the first direct determination of pt- 
broadening in semi-inclusive deep-inelastic scattering for 
charged pions and positively-charged kaons was per- 
formed on He, Ne, Kr, and Xe targets. The broadening 
was measured as a function of the atomic number A and 
the kinematic variables Q^, x or z. The broadening 
increases with A and remains constant with v, suggest- 
ing that the effect is due to the "partonic" stage and that 
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FIG. 2: From left to right, the u, Q^, x, and z dependence of (p^) for D (top row) and pt-broadening (remaining rows) for 
Tv'^ and TT^ produced on He, Ne, Kr, and Xe targets and for produced on a Xe target (bottom row). The inner error bars 
represent the statistical uncertainties; the total error bars represent the total uncertainty, evaluated as the sum in quadrature 
of statistical and systematic uncertainties. 
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color neutralization happens near the surface or outside 
the nucleus. The effect increases with increasing and 
X, and vanishes as z approaches unity. 

We thank A. Accardi, B.Z. Kopeliovich, H.J. Pirner, 
and X.N. Wang for interesting and useful discussions. 
We gratefully acknowledge the Desy management for its 
support, the staff at Desy and the collaborating institu- 
tions for their significant effort, and our national funding 
agencies for financial support. 



[1] D. de Florian, R. Sassot, and M. Stratmann, Phys. Rev. 

D 75 (2007) 114010; Phys. Rev. D 76 (2007) 074033. 
[2] L. Osborne et al, Phys. Rev. Lett. 40 (1978) 1624. 
[3] J. Ashman et al, Zeit. Phys. C 52 (1991) 1. 
[4] M. Adams et al, Phys. Rev. D 50 (1994) 1836. 
[5] A. Airapetian et al, Eur. Phys. J. C 20 (2001) 479, A. 

Airapetian et al [Hermes], Phys. Lett. B 577 (2003) 

37, A. Airapetian et al [Hermes], Phys. Rev. Lett. 96 

(2006) 162301. 
[6] A. Airapetian et al, Nucl. Phys. B 780 (2007) 1. 
[7] B.Z. Kopeliovich et al, Nucl. Phys. A 782 (2007) 224. 
[8] R. Baier et al, Nucl. Phys. B 484 (1997) 265. 



[9] R. Baier et al, Ann. Rev. Nucl. Part. Sci. 50 (2000) 37. 
[10] A. Accardi, Phys. Rev. C 76 (2007) 034902. 
[11] N.M. Agababyan et al, Phys. Atom. Nucl. 66 (2003) 
1310. 

[12] K. Ackerstaff et al, Nucl. Instr. Meth. A 417 (1998) 230. 
[13] A. Airapetian et al, Phys. Rev. D 75 (2007) 012007. 
[14] T. Sjostrand et al, Comput. Phys. Commun. 135 (2001) 
238. 

[15] B. Maiheu, Ph.D. thesis. University of Gent (2006), A. 

Hillenbrand, Ph.D thesis. University of Erlangen, DESY- 

THESIS-2005-035 (2005). 
[16] I. Akushevich, H. Bottcher, and D. Ryckbosch, 

hep-ph/9906408 
[17] R. Brun, R. Hagelberg, M. Hansroul, and J. Lassalle, 

CERN Report CERN-DD-78-2-REV (1978). 
[18] G. Ingelman, A. Edin, and J. Rathsman, Comput. Phys. 

Commun. 101 (1997) 108. 
[19] A. Airapetian et al, Phys. Rev. D 71 (2005) 012003. 
[20] A. Airapetian et al, Eur. Phys. J. C 17 (2000) 3898. 
[21] Y. Van Haarlem, Ph.D. thesis. University of Gent, 

DESY-TH ESIS-2007-033 (20 07). 
[22] A. Accardi larXiv:0808.0656k l. 

[23] B.Z. Kopeliovich, J. Nemchik, E. Predazzi, and A. 

Hayashigaki, Nucl. Phys. A 740 (2003) 211. 
[24] S. Domdey, D. Griinewald, B.Z. Kopeliovich, and H.J. 

Pirner, .arXiv:0812.283 8 [hep-ph] . 



